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i-^ \ Abstract 

Ch ' 

Q ■ In the present paper we study chemical potential of a generalized Hubbard model 

CJ I with correlated hopping at half-filling using a generalized mean-field approximation. 

For the special case of the model the approach reproduces the exact results: metal- 
insulator transition, ground state energy. Chemical potential of generelized Hubbard 
model with correlated hoppig as function of hopping integrals at zero temperature is 

\^ \ found. It is shown that chemical potential of the model is temperature-dependent. 

^? ' The dependences of chemical potential of a generalized Hubbard model with cor- 

related hopping on the model parameters are different in metallic and insulating 

0> \ phases leading to a kink at the point of metal-insulator transition in the ground 

^^ ■ state. With the increase of temperature the kink in the chemical potential curve 

jrt ■ disappears. The obtained results differ from those of the Hubbard model indicating 

the important role of correlated hopping. 

^\ 

o : 1 Model Hamiltonian 

o 

One of the simplest model describing correlation effects in narrow energy bands is the 
k> \ Hubbard model [jl| . The model Hamiltonian contains two energy parameters: the matrix 

H \ element to being the hopping integral of an electron from one site to another (to is not 

dependent on occupation of sites involved in the hopping process) and the parameter U 
being the intra-atomic Coulomb repulsion of two electrons of the opposite spins. This 
model is studied intensively (for recent reviews see Refs. 0, 0]). 

Theoretical analyses, on the one hand, and available experimental data, on the other 
hand, point out the necessity of the Hubbard model generalization by taking into ac- 
count correlated hopping. This necessity is caused by two reasons. Firstly, theoretical 
analyses [^, |^ point out the inapplicability of the Hubbard model for the description of 
real strongly correlated electron systems, in some compounds (e.g. see the estimation 
in Refs. [§, P]-[3) the matrix element of electron-electron interaction describing corre- 
lated hopping is the same order that the hopping integral or on-site Coulomb repulsion. 
Secondly, using the concept of correlated hopping and caused by it the electron-hole 
asymmetry we can interpret the peculiarities of some physical properties of narrow-band 
materials @, [ig]-[0. 



Hirsch showed that in contrast to the hopping integral of the Hubbard model (which is 
not dependent on occupation of sites involved in the hopping process) this parameter of a 
generalized Hubbard model had to depend on occupation of sites involved in the hopping 
process. Hamiltonian of the generalized in such a way Hubbard model is written as 



H = -Y.tljO'taf^ja + UY^^i^^ih (1-1) 

ijcF i 

tij = ^aa(1 - ^ia)(l - nj^) + tAB{ni^ + rij^ - 2ni^nj^) + iBBni^rij^. (1.2) 

In recent few years Hamiltonian ( |1 . 1| ) is widely used to study metal-insulator transition 



in narrow energy bands |p^-||l9| 

In Ref. [^ |T^ the necessity of the Hubbard model generalization by taking into ac- 
count the matrix element of electron-electron interaction describing intersite hoppings of 
electrons had been pointed out. The Hamiltonian of the generalized Hubbard model with 
correlated hopping is 

H = Ho + Hi + H[, (1.3) 

ifo = -/i E {Xj + Xi + 2Xf) +UY: Xl (1.4) 

i i 

H, = t{n)Y!xfXf + t{n)Y!xrxf, (1.5) 

ija ijcr 

H[ = t'{n)Y^ {v.XfX-f + h.c) , (1.6) 

ijcr 

where /i is the chemical potential, Xf' = |A;)(/| is the Hubbard operator [|20|; |0) denotes 
the state of site, which is not occupied by an electron (hole), \a) denotes the state of singly 
occupied (by an electron with spin a) i-site, |2) denotes the state of doubly occupied (by 
two electrons with the opposite spins) i-site (doublon). 

t{n)=to + nJ2Ji^kjk) =to + nTi (1.7) 

is the effective hopping integral of electrons between nearest-neighbor sites of lattice, n is 
the electron concentration, rj-^ = —1, rji = 1, 

2 

J{tkjk) = f /(/.*(r-R,)y;(r-R,)^-|v9(r'-R,)pdrdr', (1.8) 



ip{r — Rj) is the Wannier function, 

i{n) = t{n) + 2T2, t'{n) = t{n) + T2, T2 = J{iiij), (1.9) 

the prime at the sums signifies that i ^ j- 

In the model described by Hamiltonian (|1.3|) an electron hopping from one site to 
another is correlated both by the occupation of the sites involved in the hopping process 
(with the hopping integral T2) and the occupation of the nearest-neighbor sites (with 
the hopping integral Ti) which we take into account in the Hartree-Fock approximation 
(Eq. (|1.7|)). The peculiarity of model (|L3|) is the concentration dependence of the hopping 
integral t{n) in contrast to similar models. Below we shall consider the half-filling case. 



Properties of generalized Hubbard model with correlated hopping, in particular metal- 
insulator transition (MIT) has been studied in a number of recent works [17|-|ll^, pT[]-[25 
At half-filling and t' = (or t^B = 0) some exact results have been found [T^, BI 
simple cubic lattice with coordination number z MIT occurs at 



In a 



Uc = z{\tAA + tBB\) = '2z\to 



;i.io) 



If U > Uc the ground state of system is a paramagnetic Mott-Hubbard insulator with the 
concentration of doubly occupied sites d = 0, the ground state energy is equal to zero. 

For an arbitrary t' 7^ (or tAB 7^ 0) the finding of MIT criterion and the description 
of this phenomenon in a generalized Hubbard model with correlated hopping still remain 
an open problem. One of the step to solve this task is recent papers |1^, ^, ^]-|^ where 
criteria of MIT, ground state energy, concentration of doubly occupied sites have been 
found. In Refs. IT^, IT^ |2^-[E3] the authors have obtained the following criterion of MIT: 



Uc 



XIV. 



zi\tAA\ + \tBB\) = Zi\t\ + \t\) 

in agreement with the Mott's general physical ideas ^^. By means of the slave bosons 



method [27] it has been found in Ref. [^] that MIT occurs at Uc = 4:z\t + T2I; however 
here there is a problem of discrepancy of this result with the exact MIT criterion ( |1.10|) . 
The present paper is devoted to a further study of properties generalized Hubbard 
model with correlated hopping at half-filling, in particular investigation of the model 
chemical potential in the region of metal-insulator transition. 



2 Results 

The single-particle Green function in terms of the Hubbard operators reads as 

((s>p^)) = ((^r 1^;'°)) + vA{x'/\x?)) + vA{x;'\xfl) + ((xf |xf )). (2.1) 

The Green function {{Xp'^lX^P)) is given by the equation 



(i^+/i)((^ri^;')) 



■'pp' 



2t^'^^p "^ ^p^ "^ {{[Xp'^,Hi^ \Xp, )) 

+{{\x'/,H[]\x;!^)), 



(2.2) 



with [A, B] = AB — BA. Using a variant PS| of the generalized mean-field approxima- 
tion [E^ we suppose that 



yO(t tt 

^p ,-ni 



Y.<P3)X',''. X'/,H[ =Y.^,{j>3)Xf 



2 
3 ' 



(2.3) 



where e(pj) and ei(pj) are the non-operator expressions. The procedure of e{pi) and 
ei(pj) calculation is described in Ref. ^8[ (here there is a partial equivalence with the 



slave boson method |2^). Thus we obtain the closed system of equations for the Green 
functions {{Xp'^\Xp')) and {{X'^'^\Xp^)) . An analogous procedure is realized also in the 
equations for the other Green functions (p.l|). 



In this way, we find single-particle Green function and quasiparticle energy spectrum. 
For paramagnetic case in k-representation the spectrum is ||23|| : 



^i,2(k) 



-/i+ ^ T^F^, 



B{ti, - tk) -U + {4:dtiy, B = l-2d + Ad^, 

where d is the doublon concentration which is found from the equation 

+00 



(2.4) 
(2.5) 



+00 
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2iV 



E 



^k 



5k 



exp 



i?i(k) 



+ 1 



exp^ + 1 



with 



5k 



1 
2 



fi(tk-tk) U 



2Fk 



2Fk' 



1 5(tk-tk) ^ 

2 2Fk 2Fk' 



(2.6) 

(2.7) 
(2.8) 



6 = ksT, kB is the Boltzmann's constant, J\i_{E) is the spectral intensity of the Green 
function 



((^ri^p'))k = 7- 



Ak 



+ 



5k 



47r VE-Ei(k) E-E2(k) 



(2.9) 



For the Hubbard model single-particle Green function ( |2.1D and spectrum ( |2.4| ) at to = 
have the exact atomic form, and at f/ = describe the band result. For the case t' = 
quasiparticle energy spectrum ( ^.4| ) reproduces (see Ref. |2^) the exact results |T^, pi]] : 
metal-insulator transition and ground state energy. 

Chemical potential of the model is given by the equation ((Xf) = {Xf)): 



? (exp 4^ + 1 



+ 



5k 



exp^ + 1 



E 



5k 



k vexp:^^ + l 



+ 



A, 



exp^^ + 1 



(2.10) 



From Eq. (|2.1CI| ) at T = we find chemical potential of generalized Hubbard model with 
correlated hopping in the region of metal-insulator transition as 



/i 



w 



-U 



w + w 

U w — w 

2"^ 2 



{U <w + w), 
{U > w + w), 



(2.11) 
(2.12) 



with w = z\t\, w = z\t\ 



consequence of the electron-hole symmetry which is a characteristic of the model in this 
case. Note that the value oi Uc = w + w corresponds to the metal-insulator transition 
point of generalized Hubbard model. 

At arbitrary value of temperature from Eq. ( |2.1CI| ) we find the expression to calculate 
chemical potential: 



1 



1 



exp 



-E2{e) 



1 



exp -^ + 1 



de = 0, 



(2.13) 



where Ei{e), E2{e) is obtained from the respective formulae (|2.4| ) for -E'i(k), E2(k) as a 
result of the substitution ty^^e, tk^fe, t'^^^J^■ 

Fig. 1-3 where the dependences of chemical potential on the ratio U/w and tempera- 
ture are plotted show the fact that chemical potential of generalized Hubbard model with 
correlated hopping fi > U/2 and only in the atomic limit chemical potential fi = U/2. 
From Eqs. (|2.11|) and (|2.12|) one can see that chemical potential of generalized Hubbard 
model with correlated hopping depends on U, w, w. The dependences on these papram- 
eters are different in metallic and insulating phases, this leads to a kink at the point of 
MIT (see Fig. 1); t = (i.e. w = 0) and t' = 0.5t correspond to the values of correlated 
hopping parameters ri = Ti/|to| = and T2 = T2/|to| = 0.5. 

With the increase of temperature the kink in the chemical potential curve disappears 
(Fig. 2). The similar peculiarity in the evolution of free energy dependence on the pa- 
rameter U/w with temperature was noted by Mott 



From Fig. 3 one can see that in the model under consideration chemical potential 
is essentially dependent not only on the parameters w and w and also on temperature 
(in contrast to chemical potential of the Hubbard model), and what is more with the 
decrease of temperature chemical potential rapidly increases depending on the correlated 
hopping paramaters ri, T2. In high temperature region in the generalized Hubbard model 
chemical potential tends to U/2 with the increase of temperature; really, at T -h> 00 the 
probabilities of an electron finding within the lower and upper Hubbard bands are equal. 

In summary, in the present paper we have studied chemical potential of a generalized 
Hubbard model with correlated hopping at half-filling. Chemical potential of the model as 
function of the hopping integrals at zero temperature has been found. It has been shown 
that chemical potential of the model is temperature-dependent. The dependences of 
chemical potential of a generalized Hubbard model with correlated hopping on the model 
parameters are different in metallic and insulating phases leading to a kink at the point of 
metal-insulator transition in the ground state. With the increase of temperature the kink 
in the chemical potential curve disappears. The obtained result differ essentially from 
those of the Hubbard model that indicate the important role of the correlated hopping. 
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Figure 1: The dependence of chemical potential /i of generalized Hubbard model with 
correlated hopping in the ground state: the upper curve corresponds to Ti = 0, T2 = 0.5; 
the middle curve - ri = 0, T2 = 0.3; the lower curve - ti = T2 = (the Hubbard model). 
The asterisks denote the point of metal-insulator transition. 

Figure 2: Evolution of the dependence of chemical potential /i of generalized Hubbard 
model (ti = 0, T2 = 0.5) on U/w with temperature ^ = 0, 0.1, 0.3, 1 respectively. 
The lowermost curve corresponds to the behaviour of chemical potential of the Hubbard 
model. 

Figure 3: The temperature dependence of chemical potential fi of generalized Hubbard 
model with correlated hopping at U/2w = 1: the upper curve corresponds to ti = T2 = 0.3; 
the lower curve - ti = T2 = 0.1; the straight - ti = T2 = (the Hubbard model). 
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